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I. INTRODUCTION 



There has been much interests in the studies on corre- 
lated electrons in the presence of orbital degree of free- 
dom. The orbital degree of freedom is relevant to many 
transitional metal oxides It may be also relevant 

to some Ceo materials Q and samples of artificial quan- 
tum dot arrays j^]. For a theoretical understanding of 
the observed unusual properties, a SU(4) theory describ- 
ing spin systems with orbital degeneracy was proposed 
[ B|,P|. There were also numerical (lOj and perturbative 
[ llf studies of 1-dimcnsional models for these systems. 
The ground-state phase diagrams for the system with 
a symmetry breaking of SU{A) SU{2) x SU{2) were 
discussed . The phase separation was recently |l^] 

observed in experiment. Along with the rapid develop- 
ments in experiments where the metal ions has orbital 
degeneracy in addition to spin degeneracy, a theoretical 
study of such system by taking account of the kinetic 
terms caused by nearest neighbor hopping becomes im- 
portant. In a recent letter nM the Hubbard model for 
electrons with orbital degeneracy was studied. It was 
shown that the model not only has an underling SU(4) 
symmetry of spin-orbital double but also has a hidden 
charge SU(4) symmetry. An extended Lieb-Mattis trans- 
formation which maps those two SU(4) generators into 
each other is given. On the basis of elementary degener- 
ate perturbative theory, it was shown that the effective 
Hamiltonian with strong repulsive coupling at half-filling 
is equivalent to the Hamiltonian of the SO{6) Heisenberg 
model, and that at quarter-filling is equivalent to the one 
of SU(4) Heisenberg model. Some features of the model 
in one dimension was also briefly described. 

In present paper, we study the one-dimensional 
Hubbard-like model with SU(4) symmetry for electrons 
with two-fold orbital degeneracy extensively. Its exact so- 
lution is formulated by means of Bethe ansatz if sites are 
assumed to be occupied by at most two electrons. The 
features of the ground state and excited states for repul- 
sive coupling are shown. For finite N number of elec- 



trons, the configurations of quantum numbers are given 
explicitly and the spectra of excitations are obtained by 
solving the Bethe-ansatz equation numerically. For infi- 
nite N, the ground state and various kind of low-lying 
excitations are obtained on the basis of thermodynamics 
limit. 

The paper is organized as follows. In next section we 
propose the model Hamiltonian with some interpretative 
remarks. Employing standard method we carefully for- 
mulated the first quantized version of the Hamiltonian. 
We make an allowed modification so that the Bethe- 
Yang ansatz is applicable to this model. A de- 

tailed formulation from the Bethe-ansatz wave function 
to the Bethe-ansatz equation is given. In Sec. Ill, we 



explicitly show how the quantum numbers in the Bethe- 
ansatz equation should be taken for the non-degenerated 
ground state. We calculate the ground-state energy and 
Fermi momentum numerically for different numbers of 
electrons. We also compare them for different coupling 
constants. In Sec. IV, we study the excited states exten- 
sively by analyzing the possible variations in the sequence 
of quantum numbers. We indicate in each case how the 
quantum numbers change from integers to half-integers 
(or vise versa) with respect to that of the ground state. 
Numerical results of energy-momentum spectra for each 
excitation are given there. In Sec. 0, two special cases, 
weak and strong coupling are discussed. We are able 
to deduce several interesting properties from the Bethe- 



ansatz equation without solving it directly. In Sec. VI 



we consider the thermodynamics limit. After giving some 
general formulae and expressions, we study the ground 
state and calculate the g roun d-state energy explicitly for 
strong coupling. In Sec. VH, we discuss low- lying excita- 
tions in the spin-orbital sector on the basis of thermody- 
namics limit. Both contributions of holes and 2-strings 
[O are taken into account. The singlet excitatio n and 



several multiplet excitations are obtained. In Sec. VIII 



we discuss low-lying excitations in charge sector by ther- 
modynamics limit. The holon-antiholon and holon-holon 
excitations are obtained. Sec. [X is a summary of the 



main results of the paper and some conclusive discus- 
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sions. 



II. THE MODEL AND ITS EXACT SOLUTION 

We consider a Hubbard-like model for electrons with 
two-fold orbital degeneracy. The spin components are 
denoted by up ("f) and down (|), the orbital components 
by top and bottom. The four possible states of electrons 
are 



II > = 



13 > = 



T 



T 



2>= 



14 >= 



= 1 i 



i 



> 



(1) 



We use 1,2,3, and 4 to label these states from now on. 
The model Hamiltonian then reads 



(2) 

where i = 1,2,...,_L identify the lattice site, a, a' = 
1,2,3,4 specify the spin and orbital as defined in the 
above. The creates an electron with spin-orbital com- 
ponent a on site i, and riia := C^Cia is the corresponding 
number operator at site i. Eq.(y) is the Hamiltonian for 
four-component systems, and there were various discus- 
sions on multi-component Hubbard model in one dimen- 
sion Jl8|-p^. Whereas the physics that eq.(||) describes 
will be precise only when the representation space for the 
internal degree of freedom is specified ||T^ . It is specified 
to the spin and orbital, and the site is assumed to be 
occupied by at most two electrons in our present model. 

It is convenient to consider the states that span the 
Hilbert space of A^-particles 

IV' >= '/'ai,-a„(a;i,---,a;Ar)C+„^ •••C+„a«|0> ■ 

{aj},{xj} 

where xj £ {1,2,...,L}, aj G {1,2,3,4} and j — 
1,2, ...,N. The eigenvalue problem H\ip >— E\ip > be- 
comes an A'^-particle quantum mechanical problem with 
the first quantized Schrodinger operator (Hamiltonian), 

TV 

n = + f/^(l - Sa,^a,)Six,,Xj), (3) 

if site occupations of more than two electrons are ex- 
cluded where Aj^ := ■ ■ ,Xj + l, - ■ + ■ ■ ,Xj — 

The wave function of Bcthe-ansatz form in the region 
X S C{Q) :— {x\ 1 < xqi < ■ ■ ■ < xqn < L} reads 



(4) 



where x = {xi,X2, ■■■,xn), a '■= (ai, a2,'--,ajv), o,j 
stands for the spin-orbital component of the "jth" par- 
ticle; Pk is the image of a given k := (fci, k2, ■ • • ^k^) by 
a mapping P e Sn] Sn denotes the permutation group 



of A'' objects; {Pk\x) = ^jLi{Pk)j{x)j. The coefficients 
A{P, Q) are functionals on Sm^Sn- We known that any 
permutations can always be expressed as the product of 
the neighboring interchanges H^ : (• • • , Zj, Zj+i, • • •) i~> 
(• • • , Zj+i, Zj, • • •). So the requirement of antisymmetry 
for fermionic permutation is {YV'ip)a{x) = —^a{x), which 
gives 

A(P, WQ) = -7?'3J'Q0+i)^(ff p, Q), (5) 

where the spin-orbital labels are omitted and 'P'3J'Q(J+i) 
is the SU(4) spinor representation of the permutation 
operator. An immediate consequence of (H) is 
5ai,a 5{xi^ Xj)il)a{x) — 0, and hence we are allowed to con- 



sider, instead of ([5|), the following equivalent Schrodinger 
operator 



n 



N 
3 = 1 



uYS{x^,Xj), 



(6) 



here the interaction terms are independent of the spin- 
orbital labels. Then the strategy of ref. ||l^,^ can be 
used to solve the wave functions. The S-matrix that re- 
lates the coefficients A's between distinct regions in the 
configuration space of N electrons can be solved from the 
Schrodinger equation in the vicinity of hyperplane with 
{Qx)j ~ {Qx)j+i. Accordingly, we get 



sin{Pk)Qj - sin(Pfc)Q(,+i) + icVQ^'Q(J+^) 



sin(Pfc)Qj - sin(Pfc)Q(j 



+1) 



(7) 

where 2c — U/t. As it satisfies Yang-Baxter equation 
| p5[ , the Bethe-ansatz wave function is then consistently 
determined, i.e., the coefhcients ^'s in any region are de- 
termined up to an overall factor by the S^^''^'^^^^'' :— 
^'pQj.QU+i)sQO'QU+^) and that in different regions are 
related by (|). If let c — > oo in p), we know the wave 
function will be null if there are two fc's being the same 
value. So in the strong coupling limit, the fc's must take 
distinct values though there can be four states corre- 
sponding to the same kj in the absence of interaction. 

The periodic boundary condition is guaranteed pro- 
vided that A(P,7Q)e*(^'=)i^ = A{P,Q) in which 7 = 
jjTV-ijjAr-2 . . .n2ni. After applying the S-matrices suc- 
cessively, one obtain an eigenvalue equation in the SU(4) 
spinor space: 

SQi,QN gQi,Q(N~i) , , , S'^^''^^A{P, Q) 

= e-'(^^-)i^A(P,Q). (8) 

The eigenvalue problem (H) can be diagonalized by means 
of quantum inverse scattering method p5| ] by defining 
the transfer matrix as t{a) — trT{a) where T{a) — 
TAN{a - an) ■ ■ ■TA2{a - a2)TAi{a - ai), TAj{a) := 
S^^{a) e End{V^ ®V^). End means endomorphism. 
It can also be diagonalized by similar procedure as in ref. 
]l6t where the general case of continuous model with 5- 
function interaction was solved. The Bcthe-ansatz equa- 
tions reads 
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U ^-i/2(sinfcj - Xa), 



a=l 



M 



M' 



1 = -]^S_i/2(A„ - sinfci) Y[ -i(Aq - Xa') Y\_ ^-i/2{K - M/s) 

1=1 a' = l P=l 

M M' M" 

/3' = 1 7=1 



Q = l 

M' 



A/" 



/3=1 7' = 1 



(9) 



where S„(x) := [xA-inc]/[x~inc\. Eq. was given firstly 
in Ip^ . Because a particular chemical potential was in- 
troduced in the Hamiltonian, the Bethe-ansatz equation 
derived in is the same as the SU(4) Heisenberg model 
. We write out the Bethe-ansatz equation in a form so 
that it is easy to be remembered by means of the "Dynken 
diagram" [p^ of A3 Lie algebra 



the first simple root ri. The subscripts in S in eq.(|^) 
are then related to the covariant components of the 
simple roots when the simple roots are chosen as non- 
orthogonal basis, accordingly, ri — (—1/2, 1, —1/2,0), 
r2 = (0, -1/2, 1,-1/2), rg = (0,0,-1/2, 1). The high- 
est weight vector w = {wi, W2, W3) labeling the represen- 
tation of SU(4) carried out by the corresponding eigen- 
states is given by 



A 



-o 



V 

-o 



N M M' M'' 



wi^ N- 2M + M', 
W2 = M - 2M' + A/", 
= M' - 2M". 



(10) 



Where the dark dot is added to represent the charge 
rapidity kj which also has an angle of 120° relative to 



A set of coupled transcendental equations are derived 
by taking the logarithm of (||), 



M 



fcj - i ^ 0-i/2(sinfcj - Aa) 



N 



1=1 



M 



M 

81 (Aq — Aq' 


M' 


= 


-2^/„, 


a' = l 


0=1 






M' 


M" 






Qlif^P - /i/3 




- 1^7) = 


-27rJ/3, 


(3' = 1 


7=1 






M' 


M" 








Hp) + Qiii^y - 






0=1 


7' = 1 







(11) 



where 0„(x) := 2tan~-'^(^). The quantum number hj 
for charge rapidity kj takes integer or half-integer value 
depending on whether Af — 1 is odd or even. The quan- 
tum numbers /«, Jp or K-^ for flavor (we refer for the 



spin-orbital double) rapidities Aq, /i^ or !/^, take integer 
or half-integer values respectively depending on whether 
N - M - M', M - M' ~ M" or M' - M" is odd or 
even. These properties arise from the logarithm func- 



3 



tion. Once the roots are solved from the above equation 
(|ll|), the energy and momentum wiU be calculated by 

N 

E — — 2t ^ cos kj, 

N M M' M" 

1=1 Q = l 13=1 7=1 



P 



271 

T 



III. GROUND STATE 

The ground state is nondegenerate only if iV = 4n for 
n being odd numbers. This is easily seen by considering 
non-interaction case. The momentum eigenvalues (with 
periodic boundary condition) of non-interacting electrons 
are k = rn{2TT/L), m — 0,±1,±2---. For example, 
N — An for n = even, the ground state has a 70-fold 
degeneracy. In the following we will restrict ourselves to 
the case of A'o = 4n for n = odd, and consider the non- 
degenerated ground state. The non-degenerated ground 
state is a SU(4) singlet which is characterized by a 4-row 
and n-column Young tableau. The quantum numbers in 
( |rT| ) for the ground state are 



{h,} = {-{No-l)/2,... 
{14 = {-(3iVo-4)/8,. 
{J^} = {-(iVo-2)/4,.., 
{i^,}-{-(iVo-4)/8,.., 



(^0 - l)/2}, 
.,(37Vo-4)/8}, 
{No - 2)/4}, 
{No - 4)/8}. 



(13) 



Obviously, /I's and J's are consecutive half-integers while 
/'s and K's are consecutive integers. As a result of dl^ ) 
and (|l2|), the momentum of the non-degenerate ground 
state is zero. We plot the ground-state energy with re- 
spect to the filling factor N/L for various coupling con- 
stant in fig.([l]). The relation between Fermi momentum 
and the filling factor for different coupling constant is also 
plotted there. These numerical results are calculated for 
L = 20 and N from 4 to 40. For — > cx), we can take 
thermodynamics limit which will be discussed in Sec.Vl. 
The ground state for N = Nq + 1, N = Nq + 2 etc. are 
degenerate about which we will discuss later. 



IV. EXCITED STATES 

A. Excitations above the nondegenerate ground 
state 

The excited states are obtained by variation in the se- 
quence of quantum numbers {hj}, {la}, {J 13} or {K^ 
from that for the ground state. The simplest case is to 
remove one of the /I's from the sequence of ground state 
( p^ and add a new ho outside of the original sequence. 
That is 



{h,} = {- 
No -I 



No -I No -I 



no - 1, 



-I- no + 1, . 



No -I 



ho}, 



(14) 



2 ■ " ■ ' ' 2 

with \ho\ > (A^o ^ l)/2 and the other sequences in ( p^ ) 
keep unchanged. Clearly, the {No — l)/2 + no is absent in 
the set (p^. In fig.(^, we plot the numerical results for 
energy-momentum spectrum, which is a two-parameter 
family of excitation. This kind of excitations are singlet 
states. 

There are several further possibilities. After mov- 
ing one box from the fourth row in the Young tableau 
of the ground state to the first row, we get a Young 
tableau labeling a 15-dimensional irreducible represen- 
tation of SU(4) according to the knowledge of group the- 
ory. It requires M = 37Vo/4 - 1, M' = No/2 - 1 and 
M" = 7Vo/4 — 1. This causes the h's and J's to take inte- 
ger values instead of half-integer values that were taken 
for ground state. There are now M allowed values for the 
A/ - 1 distinct /'s, and M" allowed values for the M" - 1 
distinct iC's. Consequently, holes in the /'s and K's se- 
quences indispensably occurred, and then the low-energy 
states are parameterized by. 
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= K^^i + 1 - 







(15) 



where 1 < ai < 3iVo/4, 1 < 71 < No/4:. Numerical 
results of the energy-momentum spectra for this type of 
excitation are plotted in fig.(^), which is a two-parameter 
family of excitation. The states with negative momentum 
are just obtained by shifting the {hj} in (|l5| ) to the left 
by one unit. 

Moving one box from the fourth and one from the third 
row of the Young tableau for the ground state to the 
first and second row, we get the 20-fold excitation with 
M = 3iVo/4 - 1, M' = No/2 - 2 and M" = No/i - 1. 
The /'s and K's that took integer values in the ground 
state now take half-integer values; the h's that took half- 
integer values now takes integer values and the J's still 
take half-integer values. As a result, two holes in {J/3} 
are indispensably appeared, 

{h,} = {-No/2 + 1, No/2 ~ 1, iVo/2}, 
{/„} = {-3No/8 + 1, 3No/8 ~ 1}, 
{K4 = {-No/8^\,...,No/%-\}, 



Ji = - 

J/3 — Jl3 



No 



1 



+ '5l./3i. 



(16) 
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where /3 = 2,...,iVo/2 - 2 and 1 < /3i < /^a < A^o/2. 
Numerical results of the energy-momentum spectra for 
this type of excitation, a two-parameter family, are given 
in fig.(^. Similarly, the states with negative momentum 
arc just obtained by shifting the {hj} in ( |l6| ) to the left 
by one unit. 

Moving one box from the fourth and one from the third 
row of the Young tableau for the ground state to the first 
row, we get the 45-fold excitation with M ~ SNq/A — 2, 
M' = No/2 - 2 and M" = iVo/4 - 1. This makes the 
J's to take integer values instead of half-integer ones and 
K's to take half-integer values instead of integer ones. In 
this case, there exist two holes in {la} and one hole in 
{Jff}, namely 

{h,} = {-No/2 + 1,..., No/2 -I}, 
{K^} ^ {-No/8 + I,..., No/8 -I}, 

_ SNo-4 
Ji — g r Ol,ai, 

T ^0 ^ 1 i A 
Jl = ^ + 1 + t'l,/3i> 



Jf} = J;3-i + 1 + (5^,^ (/3 = 2, ^ - 2), 



(17) 



where a = 2, ...,3iVo/4 - 2, 1 < ai < as < 3iVo/4 and 
1 < Pi < No/2 — 1. Numerical results are plotted in 
fig.(|^). It is a three-parameter family of excitation. 

Taking out one box respectively from the second, the 
third and the fourth row and putting them together on 
the first row of the Young tableau, we have a 35-fold ex- 
citation. Due to M = 3No/4: - 3, M' = Nq/2 - 2 and 
M" ~ No/4:—l, we have four holes in {la}, accordingly, 

{hj} = {-No/2 + 1, No/2 - 1, iVo/2}, 
{Jp} = {-{No - 6)/4, {No - 6)/4}, 
{K^}^ {-No/8 + 1,..., No/8-1}, 

T - ^ A 

^1 — r Oi,Qi, 

4 

Ia=Io-l + l+J26a,a, (a = 2, 3iVo/4 - 3), (18) 



where 1 < l3i < < < (34 < 3A^o/4 + 1. The numer- 
ical results are plotted in fig.(^), where we did not plot 
the pattern obtained by shifting {hj} in (^8|) to the left 
by one unit, which is just the mirror image of the plotted 
pattern. This is a four-parameter family of excitation. 



M, M' and M" unchanged, the h's, J's and J's are half- 
odd integers but K^s are integers. Comparing to that for 
the non-degenerate ground state of iVo = 4n, there are 
now 37Vo/4 -I- 1 allowed values for the 37Vo/4 distinct /'s. 
So there is always a "hole" in the /'s sequence, namely. 



{hj} - {-(iVo - l)/2, (iVo - l)/2, ho}, 
/i = 37Vo/8 + <5i,ai, 

la - Ia-1 + 1 + <5a,ai (" = 2, 3iVo/4), 



(19) 



where 1 < ai < 3No/4: + 1. The J's and K's sequences 
are the same as those in (p3|) . The numerical results for 
ho > are plotted in fig.(|jy7where the zero energy corre- 
sponds to the ground state for N — No + 1. The spectra 
with negative momentum are obtained by using /ig < 0. 
It is easy to know by evaluating eq.(^0|) that each point 
in the figure represents a quadruplet. 

Adding two particles to the A^o particle ground state, 
we have N = No + 2, M = 37Vo/4 + 1, M' = No/2 and 
M" — iVo/4 for the low-energy states. This requires h's, 
J's and K's to take integer values but J's to take half of 
odd integer values. Referring to (^3|) we know that there 
must be a "hole" in the J's sequence, consequently. 



{h,} {-No/2, 
{la} = {-3No/i 
No 
4 



, No/2, ho}, 
■,3No/8}, 



Ji 



J/3 — J/ 



13-1 



1 + {(3 



,No/2), 



(20) 



No/2 + 1. The K's sequence is the 
31). The numerical results are plotted 



where 1 < Pi < 
same as that in 
in fig.(|). 

Adding three particles to the A'o ground state, we have 
N ^ No + 3, M = 3No/i + 2, M' = A^o/2 + 1 and 
M" = Ao/4 for low-energy states. This requires h's and 
K's to be half-odd integers but I's and J's to be integers, 
and then 

{hj} = {-{No + l)/2, (A^o + l)/2, ho}, 
{la} = {-{3No + 4)/8, (SA^o + 4)/8}, 
{J^} = {-Afo/8,...,iVo/4}, 
Ao 

" 8 

K^ = + 1 + S^,^, (7 = 2, A^o/4), (21) 



Ki 



Si, 



where 1 < 71 < A'o/4 + 1. The numerical results are 
plotted in fig.(||). 



B. Adding particles 

If the number of electrons are A^o + Ij Aq + 2 or A'p + 3, 
the corresponding states can be obtained by adding one, 
two or three particles into the system of Aq electrons. 
Adding one particle to the A^o ground state and leaving 



V. SPECIAL CASES 

For some special limiting cases, one is able to ob- 
tain several interesting conclusions from the Bethe-ansatz 
equation without solving it directly [^. In the follow- 
ing, we consider the weak coupling and strong coupling 
respectively. 



5 



A. Weak coupling 



Because &n{x) — + 7rsign(a;) for c ^ 0, the Bethe-ansatz 
equation (0) becomes 



kj + -Y^ sgn(sin% - X^) ^ —hj, 



a=l 



N 


M 






A/' 


^sgn(Aa - sinfc;) - 




n(Aa 




- E ^g'^(^" 


1=1 


a' = l 






/3=1 
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M' 






A/" 








- M/3') - 


f sgn(/X0 


a=l 


/3' = 1 






7=1 




M' 






Af" 








- A'/s) - 


E '^gn('^7 




/3=1 






7' = 1 



(22) 



Without loss of generality, the 7 label can be so chosen 
that is arranged in an increasing order. Then the 
fourth equation of (|2^) becomes 



M' 

E 



sgn(i^.^ - fif)) = 2is:^ + 27 - M" - 1. 



(23) 



Because \K^\ < (M" - l)/2 and M" < M' /2 (restricted 
by the Young tableau), the minimum value of right-hand 
side of (23) is — Af + 2. This means that the smallest /i/3 
is smaller than the smallest From eq.(^) we easily 
derive 



M' 

Et'^g^'^^-y- 

13=1 



-1 - M/3) - sgn(i/T, - up)] 



2(K 



7+1 



+ 1). 



(24) 



Obviously, if A'^+i — Kj — to", there must exist exactly 
to" + 1 solutions Hp satisfying < fJ^< 1^7+1- 
Again from the third equation of (P2|), we have 



M 

E 

a=l 



[sgn(^^+i - Xa) - sgn(^^ - Ac)] = 2Jp+i - 2Jp 

M" 

+2 - E['5g^i(^/3+i ~ ^7) - sgn{^ip - ^7)]. (25) 

7=1 

Obviously, if there is a such that fip < < /i/3+1 
then the right hand side of ( p5| ) equals 2(jp^i — Jp), oth- 
erwise it equals 2{Jp+i ~ Jp + 1). Proceeding the same 
strategy to the second equation and the first equation in 
( p^ ) successively, we obtain a sequence of relations that 
are summarized as follows. 



(i) If ^^7+1 — ~ m" , there exist exactly m" + 1 
solutions Hp satisfying < < ^j+i- 

(ii) For J/3+1 — Jp = to', if there exists a sat- 
isfying H0 < i^-y < M/3+I7 there will be to' A's such 
that /i^ < Aq < /i/3+1 ; otherwise if there is no v-y sat- 
isfying that, then there will be to' + 1 A's such that 

flp < Xa < /i/3+1. 

(iii) For la+i — Ja = to, if there is a /i/3 satisfying 
Aq < /i^ < Aq+1, there will exist to sin/c/'s such that 
Aq < sin/c/ < Aq+i; otherwise there will be to -f 1 sin/c/'s 
such that Aq < sinfc; < A^+i 

(iv) For hj^i — hj ~ n' , if there exists such a Aq that 
sinkj < Xa < sinfcj+i, then we will have kj^i — kj = 
27r(n' — 1)/L; otherwise we will have kj^i — kj — 2'Kn' / L. 

Apply these items to the ground state n' — m — 
to' = to" = 1, we conclude that the sequence {kj} is 
divided into groups with four successive values in each 
group. Between each pair (fc4j_3, fc4j_2), (A:4j_2, ^4j"-i) 
or (kij-ijk^j) in the same group there is one A, so to- 
tally there are three A's in each group. Furthermore there 
are two /i's each lying in between the adjacent A's. And 
between these two /x's there is always a u. 



B. Strong coupling 

For the strong coupling limit case, c ^ 00 (i.e. large U 
limit), the ration sin kj/c in the Bethe-ansatz equation is 
neglectable. The first equation of M) then becomes 



H 



2tt 

— n., - 
L 

M 

E[- 



1 



-P, 



2tan"^(2A«/c)] 



(26) 
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where rij are integers. Obviously, the Ph is just the mo- 
mentum of spin-orbital excitations in the SU(4) Heisen- 
berg chain |^,^. The other three Bethe-ansatz equa- 
tions turned to be those known for the SU(4) Heisen- 
berg chain for the scaled rapidities, Aq/c, h^/c and i>^/c. 
Equation ( ^6|) indicates that the allowed quasi-momcnta 
fc's in the strong coupling limit are quantized in units 
of 2ti/L, just like "spinless" noninteracting fermions. 
This is because the double occupancy is forbidden by 
the strong repulsive on-site coupling. The allowed quasi- 
momenta due to periodic boundary condition are deter- 
mined by spin-orbital momentum Pu- Particularly, for 
N = An with n = odd, the ground-state spin-orbital mo- 
mentum Ph is an odd multiple of tt so that the allowed 
fc's are half odd integers multiplied by 2tt/L. Therefore, 
the ground state for N — An {n = odd) is uniquely de- 
termined, i.e., it is nondegenerate. This is different from 
the noninteracting case in which the allowed fc's are al- 
ways integers multiplied by I-k jL for any N and for each 
k there are four states because of the spin-orbital degree 
of freedom. 



VI. THERMODYNAMICS LIMIT 



Replacing /c^, A^, /i^ and in eq. ( |ll[ ) by continuous 
variables fc. A, and v but keeping the summation still 
over the solution set of these roots, we can consider the 



quantum numbers /ij, /q, Jg, and as functions 
/(A), J {pi), and K{v) given by eq. (|ll|). Take /(A) as 
an example, when /(A) passes through one of the quan- 
tum numbers Jj, the corresponding A is equal to one of 
the roots Aj, similarly for J(/i), K(y), or h,(k). However, 
there may exist some integers or half-integers for which 
the corresponding A (/i, or fc) is not in the set of roots. 
Such a situation is conventionally referred as a "hole" . In 
the thermodynamics limit, — > oo, L — > oo, but N jL 
kept finite, we may introduce the density of real roots 
and the density of holes (indicated by a subscript h) , 

p(fc)+/7^fc) = {l/L)dh{k)/dk, 
C7(A)+f7,,(A) = {l/L)dI{\)/d\ 
uj{n) +uJh{^J') = il/L)dJ{n)/d^i, 
T{i^)+rhii^) = {l/L)dK{iy)/di^. 



By replacing the summations by integrals, for example 

hm -J2f(ki)= / dkp{k)f{k), 
1=1 •'^'^ 

M nB 

lim - Y^g{K)= / dAa(A).9(A), 



(27) 



and so forth, eq.(|ll 
integral equations. 



gives rise to the following coupled 



1 

p{k) + p^°\k)^ cos/c / X_i/2(sin/c- A)CT(A)dA, 

27r J^B 



(T(A)+a(°)(A) 
^(a*)+c.(°)(m) 
t{u)+t^°\v) 



Q pB pB' 

X_i/2(A-sinfc)p(A:)dfc- / Ki{X - \')(j{\')d\' - / K_^/2{X - p)Lo{p)dp, 

Q J^B J-B' 

B rB' j-B" 

K_^,^{p-X)a{\)d\- K,ip-p')Lu{p')dp' - K_^/^{p-v)T{v)dv, 



B' 
B" 



Ki{v ~ v')T{v')dv' 



(28) 



where Kn{x) := ■K~^ncl{n}(? -|- a;^), and Q, -B, B' , and 
B" in the definite integrals should be determined self- 
consistently by 



M" 

~17 



T{y)dv. 



(29) 



N 








M 


f 




1-B 


M' 


f 




J-B' 



p{k)dk, 



a{X)dX, 



uj{p)dp, 



They hold for the case in the absence of the complex 
roots. In the presence of complex roots, however, it has 
variants. In eq.(|2q) we denoted the inhomogeneous terms 



by p- 



(o) rr(°) 



and T*^"^ , which not only stand for the 



densities of holes ph, crh etc., but also the contributions 
from complex roots, two-strings. 

Once the density p{k) is solved from eq. (|2^) , we have 
the z-components of the total spin and the total orbital 
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CZ 1 

^tot _ ^ 

L ~ 2 



Q pB' 

p{k)dk + / uj{p)diJ, 
-Q J-B' 
B pB" 

a{\)dX — / T{v)dv, 



B 



B'- 
B' 



~f=2 y q'°^^^*" j_^uj{fi)df. 



(30) 



This is useful for a correct calculation of magnetizations. 
The energy is given by 



E f'^ 

— = —2t / cos kp{k)dk. 



(31) 



The highest weight vector that characterizes the corre- 
sponding representation is given by 

rQ fB rB' 

wi= L p{k)dk ~2L a{X)d\ + L (^{p.)dp, 

J-Q J-B J-B' 

/B pB pB 

(j{X)dX — 2L / uj{p)dp + L / T{i')dv, 
-B J-B' J-B" 

rB' 

(32) 



/B fB 
u!{fj.)dp — 2L T{v)dv 
-B' J-B" 



A. Ground state properties 

The ground state of the present model is a Fermi sea 
described by po{k), which is the distribution function of 
charge with respect to momentum k. The to{v) describes 
the distribution of states with spin down and orbital bot- 
tom in the i/-rapidity space. The loq{p) represents the 
distribution of either the state with spin up while or- 
bital bottom or that with spin down while orbital bot- 
tom in the /i-rapidity space. The (To(A), however, stands 
for the distribution of either state |2), |3) or |4) in the A- 
rapidity space. These distribution functions satisfy ( |2^ ) 
with B = B' = B" = oo and no holes, i.e., = 0, 
CT^"^ = 0, uj^°'> = 0, T^°^ = 0. By making Fourier trans- 
form to the second till the fourth equation of eq.(p8|), we 
have 



1 



/27r 



-c[g|/2+zg sin k 



p{k)dk 



-a{q)e 



-u){q)e 



-c\q\/2 



f(g) = a)(g)e' 



f{q)e- 



(33) 



It is not difficult to obtain a single integral equation that 
determines the po{k) 



1 cosfc „ 



sin k 



ik' 



poik')dk', 
(34) 



where kp is the Fermi momentum and 



R„{x) 



dq smh(nq) 
27r sinh(2q) ' 



,iqx-\q\/2 



Once po{k) is solved, the energy will be evaluated by the 
integral 

COS kpo{k)dk. 

-kp 

Though an explicit expression cannot be obtained from 
eq.(|3^) in the general case, it becomes easier for a nu- 
merical calculation. 

It is immediate from eq. (|33|) that the highest weight 
vector (|3^) of the ground state is a null vector. There- 
fore the ground state is a SU(4) singlet, accordingly, both 
spin and orbital are "anti-ferromagnetic" . 



B. ground-state energy for strong coupling 



The ground-state energy can be calculated explicitly 
at strong on-site coupling, c ^ 1 (i.e., U ^ t). Because 
of (sinfc — sin k') /c ^ 1 in this case and 47ri?3/2(0) = 
3 ln2+7r/2, eq. (pi) is written out up to the order 0(l/c), 



stro 1 , /ol o , TT COS_fciV 



The Fermi momentum determined from N/L 



LL Po{k)dk is 



kp = 



N 
1 



TT - (3 In 2 



TT sin(7rA^/i) 



2c 



and then the energy is calculated 
Ea _ tsm{TrN/L) 



^(^f(31n2+-) 



sin(27rA^/i) 



2t:N/L 



(35) 



where N/L is the filling factor. It becomes E^^^^'' = 
~Nt'^{6\n2 + tt)/U at half-filling N = 2L. At quarter 
filling L = N we have 



E, 



(1/4) 



.t\3 



'5 = -2— (-ln2+ -), 

N U^2 4^' 



which agrees with the result of the SU(4) Heisenberg 
model ppgd for J = 2t^ /U . Because the model is solved 
under the assumption of excluding site occupations of 
more than two, the results here are not valid for above 
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half-filling TV > 2_L in which there must exist sites oc- 
cupied by three electrons and the Bethe-ansatz wave- 
function failures at that point in configuration space. 
However, the energy is expected to be evaluate by a 
particle- hole transformation JTif , 

E{N/L, U)/L = E{A - N/L, U)/L + 3U{N/L - 2). 

(36) 

In the next section and thereafter we will study low-lying 
excitations on the basis of the thermodynamics limit. 



VII. SPIN-ORBITAL EXCITATIONS 

It is convenient to study the excitations by introduc- 
ing p{k) = po{k) + piik)/L, a{X) = r7o(A) + (JiW/L, 



uj{fj.) = llIo{p) +uji{fi)/L, and r(i/) = To(i^) +Ti{iy)/L 
where po{k), ao{X), uJo{p.), and To{h') satisfy the same 
set of integral equations as the ground state did. The 
excitation energy up to the order 0{1/L) is 



AE = - 



dk{2t cosfc + A)pi{k), 



(37) 



where A stands for the chemical potential Q can be 
replaced hy kp for a large system. Equation (|37D is valid 
for both the spin-orbital excitation and the charge exci- 
tation. The excitation energy is related to pi{k), which, 
moreover, should be solved from the following coupled 
integral equations: 



p,{k) + pi^^k) = 
a,{X) + a[''\X) = 

ri(^)+r(°)(^) = 



/OO 
i^_i/2(siii — A)cri(A)c?A, 
-OO 



K_^/2{X- smk)pi{k)dk - / Ki{X- X')ai{X')dX' ~ / K_^/2{X- ^x)uji{p)d^x 



K-i/2{p - X)ai{X)dX - / Klip - 

OO J —OO 

poo 

K-i/2{'^ - f^)'^i{f^)dp + / Ki{iy - iy')Ti{i^')diy' . 



K-i/2{i^ - v)Ti{v)diy, 



(38) 



The limits for the definite integrals are the same as that 
for ground state, which are valid for the low-lying exci- 
tations. Beyond low-lying excitations, however, the in- 
tegration limits Q, B, B' , and B" should be determined 
consistently. Here we only consider low-lying excitations. 

In order to consider the excitations above the singlet 
ground state, we must analyze the decomposition of the 
direct product of the SU(4) fundamental representation 
for N — An. Using the Young tableau, we can obtain 
that the decomposition gives rise to a direct sum of a se- 
ries of irreducible representations, i.e., (0,0,0), (1,0,1), 
(0,2,0), (2, 1,0), (4,0,0) etc. So the excitation states in 
spin-orbital sector include both the singlet (0, 0, 0) and 
the multiplets of 15-fold (1, 0, 1), of 20-fold (0, 2, 0), and 
of 45-fold (2, 1,0) or of 35-fold (4,0,0) etc. After eval- 
uating the contributions of roots and two-strings to the 
highest weight vectors that characterize the irreducible 
representations of SU(4), we can get the correct compo- 
sitions of holes and two-strings that create the possible 
excitations allowed by group theory. 



A. The multiplets 

One A hole and one i' hole together create a 15-fold 
multiplet. Let a[''\x) = cr''(A) = 6{X - A), t}°V) = 
t''(i^= (5(;/— D), and the other inhomogeneous terms in 
eq.(|3|) be null. Equation ( |38| ) is reduced to a closed form 
by Fourier transform. The excitation energy is composed 
of two terms 

Ai;(i5) = eaiX)+er{v), 

and each of them can be identified as a flavoron with 
energy 



dk{2t cos k + A)p{ (k, x), 



(39) 



where f — a, t, or lo. The Pi{k, A) is solved by 

cosfc/ (4c) 



pUk,X) 



\/2cosh[— (sinA: - A)] - 1 



cos k 



R. 



■3/2 



sin k ~ sink' 



pl{k',X)dk', (40) 



and the p\{k,v) satisfies the following integral equation: 
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pl{k,iy) 



cosfc/(4c) 



\/2cosh[ — (sinfc — ly)] + 1 
2c 



cosfc 



R. 



kp 



•3/2 



sinfc — sink' 



pl{k',P)dk'. (41) 



Two fi holes, = S{fi — fii) + 6{fi — fii) create a 

20-fold multiplet with excitation energy 

where Saj{x) is evaluated by the same integral (p9|), but 
Pi{k,fi) should solve 



cos fc/(4c) 



cosh[— (sin fc — p)] 



cos fc 



R. 



k F 



■3/2 



sin fc — smfc' 



p^ik',p)dk'. (42) 



The 45-fold multiplet is a three-hole state created by 

two A holes and one p hole, i.e., a[°\x) = S{X — Ai) -I- 

(5(A — A2) and lu[°^ (p) = 5{p — p) for which the excitation 
energy is 

A£^(45) = £a{\) + £ct(A2) -f £a;(A')- 

Four A holes create a 35-fold multiplet with excitation 
energy 

4 

^^(35) =I]£^(AJ)• 
In the above we have seen that there are three types of 
elementary excitation modes in the spin-orbital sector(let 
us call the SU(4) flavor degree of freedom). We refer these 
three elementary excitation modes the flavorons. It is 
easy to know the contributions of the holes to the high- 
est weight vectors, and to the spin and orbital. Conse- 
quently, the quadruplets (1,0,0) or (0,0, 1) are flavorons 
carrying both spin 1/2 and orbital 1/2 with energies Ca 
or Et, whereas the hexaplet (0,1,0) is a flavoron carry- 
ing either spin 1 or orbital 1 with energy e^. Clearly, 
spins and orbitals are mixed up in present isotropic on- 
site coupling. The spin orbital separation is expected 
to occur for the anisotropic cases that can be caused by 
Hund's rule. From eg. (|40|-p2|) we find that the asymp- 
totic behavior of all the three densities of roots vanish 
as the rapidities go to infinity. Thus these elementary 
excitations are gapless, i.e. e/(±c») = 0. 



B. The singlet 

By observing the contributions of two-strings to the 
highest weight vectors, we find that the flavorons can 
compound to form a singlet. In addition to placing to 



the A-rapidity a hole at A and to the z/-rapidity a hole at 
v, we take into account of three two-strings respectively 
in those three rapidities, A='= = Ao ± ic/2, p^ = po± ic/2 



and 



vq ± ic/2. In this case. 



M/L 



a{X)dX + 2, 



00 



/OO 
uj{p)dp + 2, 
-00 

/OO 
T{iy)diy + 2, 
-00 



and the inhomogeneous terms of eq.([3q) read 



p["\k)^ 



cos fcXi (sinfc — Ao), 
(A) = 6{X - A) + if3/2(A - Ao) 



+Xi/2(A_Ao)-i^i(A- Ao), 

^["Hp) = K3/2itJ- - /"o) + Ki/2{p-Po) 
-Ki(p - Ao) - Ki{p - i^o), 

t[°\v) = 5{v Kj,/2{v - vq) 

+Ki/2{v - vo) - Ki{v - pq). 



(43) 



Substituting them into eq.(38) and taking Fourier trans- 
form, we find that the term containing vq in the fourth 
equation cancels with the term containing vq in the third 
equation. After substituting the result into the sec- 
ond equation, again the terms containing po cancel each 
other. The substituting of the obtained expression of 
cri(A) into the first equation brings about an exact can- 
cellation of the terms containing Aq. As a result, the 
excitation energy is obtained 

A£;(i) ==e^(A)+£^(P), 

where Sa and e^- are evaluated by eq. (p9|) in which the 
Pi and pI satisfy the eq.(|40|) and eq.(pl|), respectively. 
Clearly, the excitations of 15- fold multiplet and the sin- 
glet are degenerate in energy. 



VIII. CHARGE EXCITATIONS 

A. The holon-antiholon excitation 

Let us consider the case of less than quarter-filling 
{N < L). we are allowed to add one "particle" out- 
side the charge Fermi sea, kp ^ [—kpjkp] but leaving 
a hole inside the charge Fermi sea, fc G [—kp,kp]. The 
calculation of the energy is required to start from 



E = -2t cos kp - 2tL 



cos kp{k)dk, 



where the integration limit kp is required to fulfill 
/Xp(fc)rffc-(A^-l)/i. 
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By introducing p{k) = po{k) + pi{k)/L etc., the exci- 
tation energy AE — E — Eq is composed of two terms 



AE{k,kp) = Bhik) + Ehikp). 
Here we introduced the holon energy 



(44) 



eh{x) ^ 2tcosx - {2t cos k + A) p1{k,x)dk, (45) 

J-kp 

and antiholon ("particle" state) energy £^(2;) — —£h{x). 
In eq.(^5|), the p1{k,x) should be solved from the follow- 
ing equation 



„,, , cosfc /sinfc — sinx 
Pi{k,x) + —R3/2 



cosfc 



R 



3/2 



c \ c 

sin k — sink' 



c 



p1{k',x)dk'. (46) 



This equation is derived from (|^ by taking 

p[''\k)=.S{k-k), 

a\°\x) = -Ki/2{X - sinfcp), 

c.i°)(^)=r(°)(.) = 0, 

which comes from the "particle" and the hole. Now we 
find that this kind of excitation consists of a holon carry- 
ing energy e{k) and an antiholon carrying energy eh{kp). 
Obviously, eq.(|4^) vanishes when k ^ kp and kp kp, 
and the holon-antiholon excitation is gapless. 



B. The holon- holon excitation 

In order to discuss states with double occupancy, we 
need to consider solutions containing complex k pairs. 
Suppose a configuration {hj} results in a complex pairs 
(two-strings), fc^ = k=F«x and two holes inside the Fermi- 
sea, fci, ^2 G [~kF,kF]- After a careful analyses of the 
Bethe-ansatz equation, one finds the string position is 
restricted to lie around a particular A solution that we 
denoted by Aq, i.e., it must satisfy 

sin(K ± ix) = Ao T ic/2 -t- ©(e"''-^). 

An exact deletion holds in the second equation of (^) for 
Ao = (sinfci+sinfc2)/2. After rewriting the Bethe-ansatz 
equation by separating the factors of the complex k pairs, 
we take the thermodynamics limit as before. In order to 
get the excitation energy we need to solve (^) with 

p^°\k)^5{k-h) + 5{k~k2) 



' cos k i4r_i/2(sm k — Aq), 



,(°) 



(A) 



{y) = 0, 



After careful calculation, we obtain the excitation energy 
AE{kuh) = Shih) + Shih) + A([/, kp). 

where Sh is the holon energy given by the same equation 
(^), and AiJI^kp) is given by 

A(t/, kp) + 2t I cos^ kKi{sin k - Xo)dk 



(2icosfc + A)pf"(fc,Ao)dfc, 



with 



pf"(fc,Ao) 
cos k '■''^ 



— kp 



COS k f sin k — Xq 
-Ri 



cos k 



R 



■3/2 



c \ c 
sin k — sink' 



pf"{k',Xo)dk' 



-ftri/2(sinfc - Ao). 



(47) 



rip)=K_,/,ip-Xo). 



Clearly the holon-holon excitation always has a gap 
Ag = 2eh{kp) + A{U,kp), which exists at any filling. 
However, the gapless modes of holon-antiholon are avail- 
able to carry charges for away from quarter-filling. It 
is easy to shown by calculating (|3^) that both holon- 
antiholon and holon-holon excitations are SU(4) singlet, 
consequently, they carry neither spins nor orbitals. 



IX. DISCUSSIONS 

In the above we have presented an extensive discus- 
sion on one-dimensional Hubbard-like model with SU(4) 
symmetry, where the sites are restricted to be occupied 
by at most two electrons.. The model was proposed to 
describe electrons with two-fold orbital degeneracy. The 
symmetries and some general features were given previ- 
ously We focused on the one dimensional case in 
this paper and studied the ground state and excitations 
by means of an exact solution. The excitation energies 
of the excited states are just sums of some particular 
terms related to quasiparticles. It provides an explic- 
itly interpretation of the separation of charge excitations 
and spin-orbital excitations. Among the charge excita- 
tions, there are gapless holon-antiholon excitations and 
holon-holon excitations with gap. Both excitations carry 
neither spins nor orbitals. They are completely decou- 
pled from the spin-orbital degree of freedom. The holons 
and antiholons move throughout the crystal at less than 
quarter-filling. Various excitations in spin-orbital sector 
consist of three basic modes which are created by the 
holes in the three rapidities for the spin-orbital double. 
That means there are three kinds of quasiparticles that 
carry spins and orbitals, i.e., two quadruplets transform- 
ing according to the fundamental or conjugate represen- 
tation, and one hexaplet forming the six-dimensional rep- 
resentation. These elementary excitations in spin-orbital 
sector are gapless. 



11 



As the on-site coupling in our model is isotropic for 
spin-orbital labels, there is no separation between spin 
and orbital. A complete separation between the spin 
waves and orbital waves is expected to occur after taking 
account of the contributions of Hund's rule. This needs 
to introduce anisotropic on-site coupling in the spin and 
orbital configuration. 

For finite N and L we plotted the excitation spectra 
by solving the Bethe-ansatz equation numerically. The 
variation of the quantum number for excited states from 
that for the ground state, and their changes from integer 
to half- integers (or vice versa) were shown in each cases. 
It provides a concrete interpretation about the collective 
excitations for the orbital degenerate electronic systems. 
The overall structure of the spectra for spin-orbital ex- 
citations changed greatly with respect to the changes of 
the correlation strength. The lowest excitation energy 
and the whole pattern are raised when the correlation 
strength decreased. However, the "particle" -hole excita- 
tion spectrum does not change much from the strong to 
the weak correlation strengths. 

In the quarter-filled band for strong repulsive on-site 
coupling, there will be no doubly occupied sites. In this 
case the total wave function will be separated into a prod- 
uct of Slater determinant of N "spinless" fermions and 
part of SU(4) Heisenberg magnets. The direct results 
from the Bethe-ansatz equation by taking strong cou- 
pling limit agree with it exactly. 

It is worthwhile to mention that the model studied 
here is not a direct SU(4) generalization of the Hubbard 
model, since a projection onto the subspace of states hav- 
ing at most two electrons at each site was made to render 
it solvable through the Bethe-ansatz. The two models are 
therefore not expected to share the same physical fea- 
tures. Considering the self-conjugate representation on 
a bipartite lattice in the strong repulsive coupling limit, 
Ref. 1^ clarified the system is dimerized with doubly 
degenerate singlet ground state and indicated the exci- 
tations are massive symmetry and antisymmetric kinks. 
In our present model, however, the local states on each 
site carry out the fundamental representation of SU(4). 

YQL acknowledges the supports of AvH Stiftung and 
interesting discussions with H Frahm. This work is sup- 
ported by NSFC- 19975040 and EYF of China Education 
Ministry. 
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FIG. 1. The relations between filling factor N/L and ground-state energy Eq (upper figure) or the Fermi momentum kp (lower 
figure). The points are calculated with U/t = 10 (□) and U/t = 2 (o) by taking N from 4 to 40 for L = 20. Noninteracting 
case is plotted by solid lines. The other lines (dashed lines for U /t = oo and dotted lines for U/t = 10) are plotted from the 
results of thermodynamics limit. 
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FIG. 2. "paxticle"-hole excitation spectrum for the model of 40 sites with 20 electrons. The overall structure of the spectrum 
does not change much between the strong (U/t = 10) and weak [U/t = 1) correlation strengths. 
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L=40, N=20, U/t=10 




L=40, N=20, U/t=2 




FIG. 3. 15-fold excitation spectrum for the model of 40 sites with 20 electrons. The overall structure of the spectrum changes 
much according to the correlation strengths. The lowest excitation energy and the whole pattern are raised when the correlation 
strength decreased. 
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L=40, N=20, U/t=10 




L=40, N=20, U/t=1 




FIG. 4. 20-fold excitation spectrum for the model of 40 sites with 20 electrons. Some of points in the upper figure are almost 
overlapped. The overall structure of the spectra changes much with respect to the correlation strength. The lowest excitation 
energy as well as the pattern are raised when the correlation strength decreased. 
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L=40, N=20, U/t=10 




FIG. 5. 45-fold excitation spectrum for the model of 40 sites with 20 electrons. 
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FIG. 6. 35-fold excitation spectrum for the model of 40 sites with 20 electrons. The complete spectra include the other part 
which is just the mirror image of the above pattern and is not plotted out. 
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FIG. 7. Excitation spectrum for one particle added into the model of 40 sites with 20 electrons. The zero energy corresponds 
to the 21 electron ground state. 
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FIG. 8. Excitation spectrum for two particles added into the model of 40 sites with 20 electrons. The zero energy corresponds 
to the 22 electron ground state. 
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FIG. 9. Excitation spectrum for three particles added into the model of 40 sites with 20 electrons. The zero energy 
sponds to the 23 electron ground state. 
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